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The recently realized multicomponent Bose-Einstein con- 
densates provide opportunities to explore the rich 
physics brought about by the spin degrees of freedom. For 
instance, we can study spin waves and phase separation Q, 
macroscopic quantum tunneling El, Rabi oscillations f|, the 
coupling between spin gradients and superfluid flow, squeezed 
spin states, vortices M and other topological excitations. 
Theoretically, there have been already some studies of the 
ground-state properties |8T~ [i"ll o f these systems and their 
line-like vortex excitations p jl2| , p^ |. In analogy with nu- 
clear physics or the quantum Hall effect, we explore here the 
possibility of observing point-like topological excitations or 
skyrmions [[l4j[il|. These are nontrivial spin textures that 
in principle can exist in a spinor Bose-Einstein condensate. 
In particular, we investigate the stability of skyrmions in a 
fictitious spin-1/2 condensate of 87 Rb atoms. We find that 
skyrmions can exist in this case only as a metastable state, 
but with a lifetime of the order of, or even longer than, the 
typical lifetime of the condensate itself. In addition to de- 
termining the size and the lifetime of the skyrmion, we also 
present its spin texture and finally briefly consider its dynam- 
ical properties. 

PACS numbers: 03.75.Fi, 03.65. Db, 05.30.Jp, 32.80.Pj 



An essential feature of a spinor Bose-Einstein conden- 
sate is that two or more hyperfine states of the atoms in 
the condensate have almost the same energy. As a result 
this spin degree of freedom becomes a relevant dynamical 
variable, which gives rise to new excitations that are not 
present in the usual single-component Bose-Einstein con- 
densates, where the spins are effectively frozen. Among 
these excitations is the skyrmion, which is a topologi- 
cal spin texture. Roughly speaking, the skyrmion is a 
point-like object that can be created out of the ground 
state, in which all the spins are aligned, by reversing the 
average spin in a finite region of space. Although topo- 
logical considerations indeed allow for these excitations, 
the fundamental question that needs to be answered is 
whether such a configuration is also energetically stable. 
A simple physical argument in favor of the stability of the 
skyrmion excitation is that, although gradients of the av- 
erage local spin increase the kinetic energy of the conden- 
sate, they also result in a reduction of the density of the 
gas and therefore in a reduction of the interaction energy. 
So we may anticipate that stability occurs when these two 
competing factors are of the same order of magnitude. 
However, when solving the appropriate Gross-Pitaevskii 
equation, we find that in principle this reduction of the 
interaction energy is always insufficient to prevent the 



skyrmion from collapsing to zero size. Fortunately, it 
turns out that for sufficiently small sizes of the skyrmion 
another stability mechanism starts to work. A number 
of atoms in the centre of the skyrmion, which we denote 
from now on as the core atoms, will be trapped by an ef- 
fective three-dimensional potential barrier, i.e., a repul- 
sive shell with a finite radius that is induced by the gra- 
dients in the spin texture of the skyrmion itself. As the 
skyrmion shrinks in size, the barrier height of the repul- 
sive shell increases and the radius decreases. This leads 
to a squeezing of the core atoms and thus to an increase 
in their energy, which ultimately stabilizes the skyrmion. 
Of course, this is not an equilibrium state of the conden- 
sate since the core atoms will tunnel over the barrier and 
give the skyrmion a finite lifetime. Before discussing this 
in detail, we first need to point out that we perform our 
calculations for a uniform condensate. In spite of this, 
our results are also valid for trapped gases, since the size 
of the skyrmion always turns out to be of the order of the 
correlation length of the condensate, which is typically 
much smaller than the size of the condensate. Moreover, 
in order to produce realistic estimates that can possibly 
be compared with future experiments, we discuss only 
the case of a spin-1/2 87 Rb condensate, since the spin-1 
23 Na condensate that has also been realized experimen- 
tally has an antiferromagnetic ground state. Considering 
other spin states or other species is in principle straight- 
forward and will be presented elsewhere. We thus con- 
sider a uniform condensate of constant density n and 
constant spinor C, z = (1,0), with all spins being oriented 
along the z-axis. This represents the ground state of 
the gas. For the skyrmion, however, both n and £ are 
position-dependent. A convenient way to take the po- 
sition dependence of £(r) into account is to express the 

spinor as £(r) = exp|— 2iQ(r) ■ SijCz, where S are the 

usual angular momentum operators for spin-1/2. The 
physical meaning of this formula is that the average spin 
at a position r is rotated by an angle 2fl (r) from its ini- 
tial orientation, which is in the positive ^-direction, with 
f2(r)/Q(r) the axis of rotation. An explicit form of f2(r) 
determines now a specific texture of the skyrmion. We 
consider only the most symmetric shape of skyrmion, be- 
cause on general grounds this is expected to have the low- 
est energy. We thus take fi(r) = u{r)r/r. The boundary 
conditions at r = and r — > oo that the function cj(r) 
should satisfy are the following. First, at r — > oo all 
spins must be oriented as in the ground state, since oth- 
erwise it requires an infinite amount of energy to create 
the skyrmion. This implies lim r _ >00 o;(r) = 0. Along the 
z-axis the spins are also not rotated by our ansatz for 
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fi(r), so in order to have a nonsingular texture of the 
spinor with a nonzero winding number, we must require 
that w(0) = 2n. Finally, to avoid a singular behavior of 
fi(r) itself, we take only functions w(r) with zero slope 
at the origin. In summary, to(r) is therefore a monotoni- 
cally decreasing function that starts from 2tt at the origin 
and reaches zero when r — > oo. The specific functional 
form of uj(r) at intermediate distances between zero and 
infinity is not crucial for the stability of the skyrmion. 
Only the above boundary conditions are important for 
that. From a quantum mechanical point of view, the 
condensate is described by a macroscopic wave function, 
or order parameter, ip(r) = ^/n(r)£(r). Furthermore, 
the grand-canonical energy of the gas can in the usual 
mean-field approximation be expressed as a functional of 
n(r) and £(r) M. In detail we have 



£?[n(r),C(r)] = / d(r 
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where m is the mass of the atoms, \i is their chemical po- 
tential and T 2B = 4nh 2 /2m is the appropriate coupling 
constant that represents the strength of the inter-atomic 
interactions in terms of the positive scattering length a. 
It is clear from the above expression that the gradients 
in the spin texture lead to an energy contribution, which 
is proportional to |VC(r)| 2 . Using the above-mentioned 
form of C(r) and inserting the explicit forms of the spin- 
1/2 matrices, the square of the spinor gradient can be 
written explicitly as 



|VC(r) 



Kr)) 



dr 
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In principle, both n(r) and t-u(r) can be calculated ex- 
actly by minimizing the energy functional in Eq. (|^) 
with respect to arbitrary functions n(r) and u>(r) that 
satisfy the boundary conditions. However, the result- 
ing nonlinear and coupled equations are quite difficult 
to solve. Therefore, we use here a variational approach 
and take for the simple ansatz uj(r) — 4cot~ 2 [(r/A) 2 ] , 
where the variational parameter A physically corresponds 
to the size of the skyrmion. This ansatz is chosen here 
because it automatically incorporates the correct bound- 
ary conditions and leads to a minimal skyrmion energy 
as compared to various other ansatzes that we have con- 
sidered. Having specified a->(r), we then calculate n(r) 
exactly by solving numerically the differential equation 
for n(r) obtained by varying E[n(r), £(r)] with respect 
to n(r). Substituting this density profile back into the 
energy functional, the energy of the skyrmion becomes 
a function of A only and the equilibrium properties can 
be obtained by minimizing this energy with respect to A. 
Inserting our ansatz for uj(r) in Eq. (§), ft 2 |VC(r)| 2 /2m 
takes the shape of an off-centered potential barrier with a 
maximum of 24.3ft 2 /2mA 2 located at 0.68A. We now dis- 
tinguish between two cases. The first case occurs when 



the height of the barrier is lower than the chemical poten- 
tial /i of the atoms. In this case, atoms can move freely 
across the barrier. In the second case, the barrier height 
is higher than the chemical potential and thus atoms be- 
come trapped behind the barrier near the centre of the 
skyrmion. Equating ?i 2 |V£(r)| 2 /2m and fi gives therefore 
the maximum value of X ma x ~ 5£ below which trapping 
takes place. Here £ is the correlation length given by 
£ = 1/ VSnan, and n is the density at r 3> A. It should be 
noted that to obtain this value for X max , we used the fact 
that at sufficiently large radial distances, gradient terms 
in the equation of motion derived from the energy func- 



tional vanish, and wc obtain fj, = T 
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For A < A,, 



we calculate first the equilibrium size of the skyrmion by 
minimizing for a fixed number of core atoms the total en- 
ergy of the condensate. Next we also calculate the tunnel- 
ing rate for the core atoms to escape to the outer region 
(see Appendix) . The results of these calculations are pre- 
sented in Fig. 1, where we plot the equilibrium size of the 
skyrmion as a function of the number of the core atoms 
trapped by the texture barrier, together with the corre- 
sponding tunneling rate. As mentioned previously, we 
use the parameters of 87 Rb with a density of 10 13 cm" 3 . 
We see from this figure that only a few atoms in the core 
of the skyrmion are needed to stabilize it and to give it a 
sufficiently long lifetime. We note that the tunneling rate 
is considerably less than the decay rate of the condensate, 
which is due to two-body collisions and equal to Jlj| Gn 
with a measur ed value of § G « 2.2 x 10 14 cm 3 /s. 
We also mention that the lifetime becomes even much 
larger for slightly smaller values of n. Finally, we present 
some of the interesting properties of the skyrmion, which 
are its texture and its dynamics. The texture can be 
best presented in terms of the three average spin com- 
ponents (S x ) (r) = tC(r)&C(r), (S y ) (r), and (S z ) (r) in 
the three Cartesian planes. These are shown in Fig. 2. 
The most interesting dynamical property of the skyrmion 
comes from the fact that if all spins of the texture are 
rotated around the z-axis by a constant angle 9, i.e., 
£(r) — ^ exp (— i9S z )C( r )i the resulting skyrmion will have 
the same energy. As a result the angle 9 undergoes phase 
diffusion. It also leads to Joscphson-likc coupling be- 
tween two skyrmions, which will have important conse- 
quences for the physics of a skyrmion lattice An- 
other dynamical property is the center of mass motion 
of the skyrmion, which can be shown to be identical to 
that of a free particle. Both dynamical properties will be 
discussed in detail in a future publication (U. and H. T. 
C, in preparation). 



APPENDIX 

To calculate the energy of the skyrmion we need to 
solve the equation for the density profile that is obtained 
from minimizing the energy functional in Eq. (||). We 
solve this equation numerically for the region outside the 
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core. Inside the core we solve the equation analytically by 
using the Thomas-Fermi approximation, which amounts 
to neglecting the gradients of the density profile. Using 
our ansatz for u>(r), the texture gradient potential V(r) — 
S 2 |VC(r)| 2 /2m reads 



V(r) 



32ft 2 (r/A) 2 [3 + 2(r/A) 4 + 3(r/A) s 
2mA 2 [i + (r/A) 4 ] 4 



(3) 



For small r/A this potential can be approximated by a 
harmonic potential with a characteristic frequency loo = 
•\/96ft 2 /2m 2 A 4 and width I = A/%/96, as shown in the 
dotted curve in Fig. 3. The use of a Thomas- Fermi ap- 
proximation is justified when the mean-field interaction 
energy is larger than the spacing between the lowest en- 
ergy levels of the harmonic trap. Specifically, the ratio 
2Na/l = 2y/96Na/X should be bigger than 1. From Fig. 
1, we observe that this ratio equals approximately 1 for 
N = 4 and it increases for larger N. The lifetime of the 
skyrmion is estimated by calculating the tunneling rate 
from the core to the outer region over the barrier V(r). 
To this end we employ the following WKB expression for 
the tunneling rate jl8| 



2tt 



exp 



-2 / dr 
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where [i C ore is the chemical potential of the core atoms. 
The radial points T\ and Ti are the points where v(r) 
and fJ, cor e intersect as shown in Fig. 3. The chemical 
potential /i core is calculated by differentiating the total 
energy of the core with respect to the number of core 
atoms. 
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FIGURE CAPTIONS 



FIG. 1. Decay rate and size of the skyrmion as a function 
of the number of atoms trapped in the core. The two curves 
are calculated for the parameters of a 87 Rb spinor condensate 
with a density of 10 13 cm -3 . The size is calculated in units 
of the correlation length , which for these parameters equals 
0.85^im. 



FIG. 2. The average spin texture of the skyrmion. Shown 
are the components (S z (r)) and (S y (r)) in the three Cartesian 
planes. A value of A ~ £, which corresponds to 20 core atoms, 
was chosen to compute these quantities. The distances are 
in units of the coherence length £. The (S x (r)) component 
can be obtained from the {S y (r)) component, by using the 
cylindrical symmetry of the skyrmion. Fig. la shows most 
clearly that our skyrmion is composed of two coaxial tori. The 
two pink rings represent the cores of these two tori where 
(S z (r)) reaches its maximum negative value -1/2. Fig. lb 
shows the cross sections of the two tori. 



FIG. 3. Figure 3 The potential barrier produced by the 
skyrmion texture. The value of A used for this figure is again 
approximately xi, which corresponds to 20 core atoms. The 
lifetime of the skyrmion is determined by the tunneling of 
core atoms with a chemical potential fj, CO re to the right of the 
barrier. The square root of the shaded area is the integrand 
of Eq. ^ . The dotted curve represents the harmonic approx- 
imation to V(r) for small r/A. The inset shows u(r) for the 
same value of A. 
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